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We discuss in detail the shear viscosity coefficient r\ and the viscosity to entropy density ratio n/s 
of a hadronic gas comprised of pions and nucleons. In particular, we study the effects of baryon 
chemical potential on n and n/s. We solve the relativistic quantum Boltzmann equations with binary 
collisions (ivn, nN, and NN) for a state slightly deviated from thermal equilibrium at temperature 
T and baryon chemical potential fj,. The use of phenomenological amplitudes in the collision terms, 
which are constructed to reproduce experimental data, greatly helps to extend the validity region 
in the T-fi plane. The total viscosity coefficient n(T, fi) = n^ + n N increases as a function of T and 
fi, indirectly reflecting energy dependences of binary cross sections. The increase in /j, direction is 
due to enhancement of the nucleon contribution n N while the pion contribution rf diminishes with 
increasing [i. On the other hand, due to rapid growth of entropy density, the ratio n/s becomes 
a decreasing function of T and /x in a wide region of the T-fJ. plane. In the kinematical region we 
investigated T < 180 MeV, n < 1 GeV, the smallest value of n/s is about 0.3. Thus, it never violates 
the conjectured lower bound n/s = l/4n ~ 0.1. The smallness of n/s in the hadronic phase and its 
continuity at T ~ T c (at least for crossover at small fi) implies that the ratio will be small enough 
in the deconfined phase T <; T c . There is a nontrivial structure at low temperature and at around 
normal nuclear density. We examine its possible interpretation as the liquid-gas phase transition. 



INTRODUCTION 



Shear viscosity of a hot QCD matter has been attract- 
ing much attention in recent years. The major reason 
for that is the intriguing experimental discovery that 
the matter created in heavy-ion collisions at Relativis- 
tic Heavy Ion Collider (RHIC) in Brookhaven National 
Laboratory could be close to a perfect fluid [l[ . This un- 
expected result has driven people to think about strongly- 
interacting quark-gluon plasma (abbreviated as "sQGP") 
which may be realized at temperature just above the crit- 
ical temperature T <^ T c . As a result of strong-coupling 
nature, sQGP is thought to have very small shear viscos- 
ity, which is however not directly confirmed yet in a sat- 
isfactory way. In fact, only a few things are understood 
about sQGP because we will not be able to investigate 
it within standard perturbative QCD techniques. The 
only technique available now (except for lattice simula- 
tions which are not analytic methods) is the one based 
on the AdS/CFT correspondence, which relates strongly- 
coupled supersymmetric Yang-Mills theories to weakly- 
coupled gravity theories. 

There is an interesting outcome from the AdS/CFT 
analysis in relation to the shear viscosity: It has been con- 
jectured that there would be a lower bound in the "shear 
viscosity coefficient to the entropy density ratio" (or sim- 
ply the "viscosity to entropy ratio") n/s > l/4ir Q. We 
call this "the KSS bound" after authors' names of Ref. 0]. 
The lowest value n/s = 1/A-k is satisfied by several super 
Yang-Mills theories in the large N c limit (strong coupling 
limit), which suggests that the bound could be universal. 
Of course, there is no guarantee for this bound to hold 
in real (non-supersymmetric) QCD whose gravity dual 



is not found, but interestingly enough, the values of n/s 
extracted from RHIC experiments [3| and from lattice 
simulations [J, [5( seem to be small enough and close to 
the lower bound. 

On the other hand, there is an important empirical ob- 
servation which can be seen in many substances such as 
helium, nitrogen, and water: The ratio n/s has a min- 
imum at or near the critical temperature [|| (see also 
[7|]). More precisely, the ratio shows a cusp at T c for the 
first order transition, while it has a convex shape for the 
crossover with its bottom around the (pseudo) critical 
temperature. Since this behavior is observed in many 
substances, it is expected to be universal. Recall that 
the phase transition in QCD is most probably crossover 
at least for low densities. Therefore, what we naturally 
expect for the QCD matter from the two observations 
mentioned above is that the ratio n/s will have the min- 
imum at T ~ T c , and the numerical value at that point 
will be close to the KSS bound n/s ~ 0.1. It is of primary 
importance to check whether this expectation is indeed 
the case or not, and to understand the properties of the 
QCD matter around T c not only for T ^T c . These con- 
siderations motivated us to investigate the shear viscosity 
in QCD from the hadronic phase T ^ T c . Notice that we 
can indirectly study the properties of sQGP from below 
T c because physical quantities such as the ratio n/s will 
be continuous at T c for the crossover transition. 

Transport properties of a meson gas have been stud- 
ied by several people. Many of the calculations arc 
based on the Boltzmann equations with the Chapman- 
Enskog method which is a standard approach for weak 
dissipative phenomena, especially for computing trans- 
port coefficients [§]. Differences among several papers 
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@, 03 EH d, include kinetic or statistical proper- 
ties of particles (relativistic or nonrelativistic, quantum 
or classical), species of mesons (pions, kaons, etc.), and 
the cross sections in the collision terms. For example, 
Ref. [l2l ] treated a nonrelativistic quantum pion gas with 
the binary cross section given by the leading order chiral 
perturbation theory (LO-ChPT). Recently, similar prob- 
lems have been revisited in relation to the KSS bound 
[HI, [Hj]. In Ref. [lH], the ratio 77/s computed with the 
cross section in LO-ChPT turned out to violate the KSS 
bound for temperature beyond T c ~ 170 MeV, and it was 
speculated that such violation could be related to the ex- 
istence of phase transition. However, soon after that, it 
was shown in Ref. [HI that the KSS bound is not violated 
in a pionic gas if one computes the shear viscosity with 
a phenomenological cross section using the experimental 
phase shifts. What we have learned from these papers 
is the following: (i) we have to be careful when we use 
the cross section from effective field theories, and (ii) the 
ratio 77/s of a relativistic pion gas is small enough at rel- 
atively large temperature T ^ T c , but does not violate 
the KSS bound 77/s > 1/4tt. 

At this point, there comes a natural question: how 
does 77/s change if one adds nucleons to the pure pion 
gas? Naively, we expect that the ratio 77/s will decrease 
as number of nucleons is increased because the pion cross 
section will effectively enhance in the presence of nucle- 
ons, yielding smaller shear viscosity, while the entropy 
will increase. Thus, there is a chance that the ratio could 
violate the KSS bound if the bound does not change. No- 
tice that the pion-nucleon gas is the minimum requisite 
which allows us to study the effects of baryon chemical 
potential /i. Therefore, it is quite interesting and impor- 
tant to investigate the fi dependence of 77/s in the pion- 
nucleon gas. Such investigation will also urge people to 
study the (possible) /_t-dependence of the KSS bound. In 
fact, the "universality" of the KSS bound has not been 
tested at finite baryon chemical potential. 

Most recently, the authors of Ref. [14[ have applied 
their framework to the pion-nucleon gas to study the be- 
havior of 77/s in the T-\i plane [lijj . However, their fo- 
cus was not on the KSS bound but on the new finding: 
a valley structure in 77/s at low temperature and large 
chemical potential which they argued would correspond 
to the nuclear liquid-gas phase transition. Although this 
is a very interesting suggestion, their results should be 
critically checked since they are based on the effective 
field theories whose validity region is severely limited. 
On the other hand, a realistic calculation of the viscos- 
ity in the pion-nucleon gas was performed some time ago 
by Prakash et al. [l(|. Remarkably, they used the bi- 
nary cross sections in the collision terms, which roughly 
reproduce experimental data. However, unfortunately, 
the dependence on baryon chemical potential was not in- 
vestigated in detail. Besides, this calculation is based 
on classical Boltzmann equations, and thus cannot be 
applied to relatively large chemical potential where the 
effects of Fermi statistics is expected to be large. 

In view of the present situation mentioned above, what 



we should do is rather evident: for the purpose of study- 
ing the fj, dependence of the shear viscosity 77 and the 
ratio 77/s in a pion-nucleon gas, we treat the relativistic 
quantum Boltzmann equations with binary cross sections 
which are determined to reproduce experimental data. 
We are very careful about the range of validity of our 
framework. We also check whether the valley structure 
found in Ref. [16| | indeed exists even with the phenomeno- 
logical cross sections. It is also important to compare our 
results with those from hadron cascade simulations. For 
example, the shear viscosity coefficient is computed for 
a meson gas in Ref. ljj and for a meson-baryon gas in 



Ref. [18(. 

The paper is organized as follows: in the next section, 
we explain the relativistic quantum Boltzmann equations 
for a dilute pion-nucleon gas, and define the cross sec- 
tions we use in the collision terms. We treat the small 
deviation from the thermal equilibrium to the linear or- 
der (the Chapman- Enskog method), and give the shear 
viscosity coefficient 77 through the solutions to the Boltz- 
mann equations. In Sect. Ill, we present our numerical 
results for a pure pion gas and a pion-nucleon gas mix- 
ture. We introduce a criterion which measures the valid- 
ity region of the calculations. The use of phenomenolog- 
ical cross sections is very important to enlarge the range 
of validity. We discuss in detail the effects of chemical 
potential on 77 and 77/s and examine the interpretation 
of the valley structure as the liquid-gas phase transition. 
Summary is given in the last section. 



II. KINETIC THEORY OF A HADRONIC GAS 
MIXTURE 

A. Quantum Boltzmann equations 

We first explain our theoretical framework which is 
necessary for computing the shear viscosity coefficient. 
Consider a dilute gas of pions (it) and nucleons (N) in 
which particles interact with each other through binary 
collisions. Noncquilibrium processes such as relaxation 
to thermal equilibrium can be described by kinetic equa- 
tions for one particle distribution functions /^(x, p,t) 
and f N (x,p,t) (below, we suppress x-dependence) . For 
simplicity, we assume that the gas is isospin symmetric, 
and thus /"(p, t) and f N (p, t) are isospin averaged distri- 
butions. The relativistic quantum Boltzmann equations 
(more precisely, the Uehling-Uhlenbeck equations) of this 
hadronic gas mixture are then given by 



-a M r(p) = c™ lt, n + c« N [r, f N ] , (i) 



j^d»f N (p)=c NN [f N j N ] +c N * [r,r] , (2) 



where E£< N = 
defined as 



7T}iV 



p 2 and the collision terms are 
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c n* +c nN = 



^ J dv™[fU2 (i + /? ) (i + /;) - (i + /n (i + fi) /?/;} 

5JV | dr^{/f/j (i - f 3 N ) (i + /;) - (i - /f) (i + /f) /f /;} , 
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(4) 



We have used shorthand notation f*' N = 
f n ' N (ki), f£> N = f^ N (p) and g^, g N are the de- 
generacy factors g n — 3, gN = 2. For the collisions 
between the same species (tttt, NN), we have added a 
factor 1/2. The factors (1 + f 71 ) and (1 - f N ) represent 
the Bose-Einstein and Fermi statistics of particles, 
respectively. Finally, dT in the integrants are invariant 
measures with the scattering amplitudes squared: For 
example, 



for 7T and N, 



dT 



TTiV 



N 



(2TT) 4 S^(k 1 +k 2 
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P) -Q d3k i 



(2tt) 3 ' 

where M^n is the elastic scattering amplitude for 
N(ki) + 7r(fc 2 ) -> N(k 3 ) + ir(p). Explicit form of the 
scattering amplitudes will be shown at the end of this 
section. These are the basic ingredients of the kinetic 
theory. 

Before we discuss how to solve the Boltzmann equa- 
tions ([1]) and ((2]), we need to know the equilibrium states. 
We define them without solving the full Boltzmann equa- 
tions: They are given by the distributions f£, which 
make the collision terms vanish. Namely, C 7r7r [/o r , ffi] + 
C^lfSJ?) = C N "lf N JS}+C NN {f N J N } = 0. These 
conditions are easily satisfied by the following Bose- 
Einstein and Fermi-Dirac distributions if the common 
temperature T — 1/(3 and hydrodynamic velocity V M 
are used 



fo N ( P ) 



1 



e 0V^ _ i ' 
1 

e /3(V t ,p^-M) -|- 1 



(G) 
(7) 



where [i is the baryon chemical potential. The parame- 
ters T, V 11 and /j, can, in principle, depend on the coor- 
dinates (local equilibrium), but when we compute quan- 
tities in thermal equilibrium, we simply select the rest 
frame V = (1, 0, 0, 0) so that V^p* = E. 

In evaluating the entropy density, we use the expres- 
sion in the equilibrium state as is done in the literature 
because the deviation from the equilibrium is assumed to 
be small. Namely, by using the grand partition functions 



lnZ„. = -Vg v 



d 3 p 
(2^)3 



In 1 - e T 



({' n , )■■ , 

lnZ N = Vg N j ^yj ln l 1 <-"^~ 



(8) 
(9) 



one obtains the total entropy density s = s w + Sjv a s 
follows: 

1 d , 

S-K = —-^TYVlZt; 
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B. Shear viscosity coefficient 

The shear viscosity coefficient is defined through the 
deviation of spatial components of the energy momentum 
tensor in the linear order with respect to fluctuation from 
the equilibrium. Consider a nonequilibrium state which 
is slightly deviated from the global equilibrium. Small 
deviation of the space components of energy momentum 
tensor (T ij 
and traceless parts 



T£> + ST ij ) can be divided into traceful 



61 'i = c;rt"'V-Vi 2;/{V'l 

i / s , * a <y J 



hi ' 



(10) 



(v l w) trl = - (vV' + v'v* -— v-v, (ii) 



where V 1 is a space component of the hydrodynamic four 
velocity (which is common for 7r and N) 

jd^ N {p) ' 

Eq. (I10|) is the definition of the shear and bulk viscosity 
coefficients r\ and £. The flow vector V 1 , as well as T 
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and /i, is in principle arbitrary and can depend on spa- 
tial coordinates. But below we consider the case where 
only the flow vector V 1 depends on the coordinates, and 
in particular, its divergence is vanishing: V l V l = and 
\j l yi zfz (i ^ j). This situation corresponds to consid- 
ering only the shear viscosity coefficient rj and ignoring 
all the other transport coefficients such as bulk viscosity 
or heat conductivity. 

If one knows distributions f^ :N (p) = fo' (p) + 
&f*' N {p) as the solutions to the Boltzmann equations 
([1]) and @, one can explicitly compute the (deviation 
of) energy momentum tensor: 



ST 3 



d3p plp3 -5f N ( P ). 



(2tt)3 E? 



(12) 



As mentioned above, we consider only the deviation 
fifn.N originates from the shear. Then, it is quite 
convenient to parametrize <5/ 7r and Sf N as follows (pi = 
Pi/p)- 

6r = -fZ(l + f5)0B*(p) 



Sf 



N _ 



f»(l-f»)/3B»(p) 



(14) 



where we have introduced new quantities B*' (p) to be 
determined by the Boltzmann equations. 

Substituting Eqs. {T3J) and QU) into Eq. (TS]) and com- 
paring the result with Eq. (p~0|) . one finds the shear vis- 
cosity coefficient rj as a function of unknown functions 
B n > N (p): 



V 



9* (3 
15 

9n(3 
15 



d 3 P KV + fS)* 



E n 
^p 



(2tt) 3 
d 3 P f N (l 



p l B%p) 



^P 



p 2 B N {p). (15) 



The unknown functions B 7r ' N (p) (or equivalently, the de- 
viations 6f*' N ) are numerically determined by solving 
the Boltzmann equations that are linearized with respect 
to df*' N . This procedure corresponds to the lowest order 
Chapman-Enskog method. Here we discuss only the out- 
line of the procedure to solve the Boltzmann equations. 
More details are discussed in Appendix B. 

After the linearization, the Boltzmann equations JT]), 
([2]) become a coupled linear equations for B 7r (p) and 
B N (p). Following Ref. ijj], we solve these equations in 
the functional space spanned by the orthogonal polyno- 
mials {W(„)(p), n = 0,1,2,...}. Let us expand B rr ' N (p) 
by these bases: 



B*> N (p) 



n=0 



(16) 



where WT\ (p) is a polynomial of the order n, and bj^ is 
the coefficient independent of p. Notice that W^(p) and 
W^(p) are not equivalent to each other. Indeed, they 
are defined so that they satisfy the following different 
orthogonal conditions: 

d 3 P /oXl+ZoO 2 TT/7T / UI/7T / \ £ tit 

p W {n) (p)W [m) (p) = 5 nm L (n) , 



(2tt 



^p 



(17) 



d 3 P fo N (l-f N ) 
» 3 



^P 



p'W^ ) {p)W^ ) {p)=8 nm L^ ) 



(18) 

The normalization factors LJ n s and are not chosen 
to be 1. Instead, we choose the polynomial so that the 
coefficient of the term with the highest degree is unity 
(such polynomials are called "monic"). Namely, the first 
three polynomials have the following form: 

W(„)(p) = 1, 
W(i)(p) = P + ci, 
W(2)(p) = p 2 + d lP + d 2 . 

Parameters c±, d\ and d 2 are uniquely determined by 
the orthogonal conditions (thus independent of the dy- 
namics). In fact, Wi n ) nas n unknown parameters, which 
are uniquely determined by n orthogonal conditions with 
lower polynomials (W( m ), m = 0, • •• , n— 1). Thus, LJ n ^ 

and in Eqs. (fT7|) and (fT5)) are known after we com- 
pletely determine Wi n y Practically, the expansion §TE\i 
is well approximated by the first few terms. Thus, in the 
present analysis, we take only the first three terms: 



fl*(p) : 
B N ( P ) 



h + b h w (D^ + b h w (2)(p^ 

+ bf 1) W { N 1) (p) + b^ ) W^(p), 



h N 



(19) 
(20) 



where we have used the definition WZ,\ N = 1. The coef- 



ficients bJ\ are numerically determined. Once we know 
these coefficients, we can compute the shear viscosity co- 
efficient rj from Eq. (fl~5|) . We have checked that the re- 
sults do not change even if we take up to the fourth terms 
(n < 3). 

Two comments are in order about our formulation. 
Firstly, we recall that the expansion of B 7r ' N (p) in 
Eq. p6p is a familiar technique in solving the Boltzmann 
equation by the Chapman-Enskog method. If one treated 
a classical Boltzmann equation, measure of the orthogo- 
nal condition would be given by a much simpler distribu- 
tion, the Maxwell-Boltzmann distribution. In this case, 
the polynomials that satisfy the orthogonal condition are 
given by famous functions, the Sonine polynomials (l9j |. 
For the quantum Boltzmann equation, however, the mea- 
sure in the orthogonal condition is given by either the 
Fermi-Dirac or Bose-Einstein distribution, as shown in 
Eqs. (|17p and (fT5)) . In this case, the polynomials satisfy- 
ing them are not known, and we have to find them order 
by order. 
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Secondly, note that the shear viscosity coefficient r\ de- 
fined by Eq. (I15|) may be formally written as 



r) = r] + r) 



N 



(21) 



It apparently looks 'additive' with respect to each con- 
tribution. Indeed, if one considers a pion gas system 
without nucleons, one finds exactly the same expression 
as the first term of Eq. (fTS)) . One might then be tempted 
to conclude from the expression (|2"Tj) that the inclusion 
of nucleons always contributes to enhance the value of 
shear viscosity. However, such argument does not make 
sense because the function B 7r (p) itself will change by the 
inclusion of nucleons. We will see later that rf indeed de- 
creases as the effects of nucleons become large (i.e., with 
increasing chemical potential). 



C. Scattering amplitudes of binary collisions 

Let us show the explicit expression of the scattering 
amplitudes in the collision terms ©, Q, and (5J). Note 
that we have assumed factorization of the scattering am- 
plitudes from the products of one-particle distributions in 
the collision terms. This is physically natural in a dilute 
gas where each collision is simply treated as an indepen- 
dent binary collision (ignoring higher order multiparticle 
correlations). Therefore, as long as we consider a dilute 
gas system where the Boltzmann equation is applicable, 
it is reasonable to use the amplitudes for two particle 
scatterings in the vacuum (free space). This is a great 
merit in computing transport coefficients. If one follows 
the microscopic Kubo formula to compute the transport 
coefficients, it is quite nontrivial how to include the ef- 
fects of physical processes. On the other hand, if one 
uses the Boltzmann equation, one can easily incorporate 
the physical cross sections in the collision terms, which 
is however less justified from the first principle. In the 
present paper, we examine two different parametrizations 
for the scattering amplitudes. One is based on theoretical 
calculations, while the other is constructed from the ex- 
perimental data. More precisely, we use the amplitudes 
from the low energy effective field theory (EFT) on the 
one hand, and the phenomenological amplitudes designed 
to reproduce experimental data of the elastic scatterings 
on the other hand. 

Scattering amplitudes based on EFT are the following. 
First of all, the isospin-averaged tttt scattering amplitude 
is given by the leading order Chiral Perturbation Theory 
(LO-ChPT) 



9/ 4 



21mt + 9s 2 -24m 2 s + 3(t 



u) 2 }, (22) 



where / ff is the pion decay constant (/V = 93 MeV), s, t 
and u are the Mandelstam variables for the scattering 
Tr(ki) + 7r(/c 2 ) — > 7r(fc 3 ) + 7r(p). 

Next, for the irN scattering, we use the results of low 
energy effective theory (LO heavy baryon ChPT) [211] . 



The isospin averaged nN scattering amplitude in the 
Center-of-Mass (CM) frame is 



|M T 



N 



(2m N ) 



! {4| 



5-1 



9cm \ h+ 



| 2 }, (23) 



where q cm and 9 cm are, respectively, the magnitude of 
pion momentum, and the scattering angle in the CM 
frame. Two functions g- and h + arc 



9- 



h+ = 



4t- ( 2 - 2 



2m 2 



2/ 2 



where Qa = 1-26 is the nucleon axial charge, and u> is 
the pion energy in the CM frame. Note that the overall 
factor in Eq. (|2"3"|) is not the same as in Ref. [l|| though 
both are based on Ref. |2lj . This is because we have 
used the standard normalization factor for the spinors 
u r (p)u s (p) = 2mNS rs instead of u r (p)u s (p) = {(E N + 
mjv) /2rriN}S rs which was adopted in Ref. [2l| . 

Lastly, let us consider the NN scattering. In the low 
energy EFT, scattering amplitudes in the CM frame for 
fixed spin and isospin can be parametrized in terms of 
the scattering length a and the effective range r as 



| M N N , | 2 = 64tt 2 s- 



1 



1 2 

2 ™cm 



iq c 



(24) 



where q CI[L is the magnitude of nucleon momentum in the 
CM frame. For each process, parameters a and r are 
determined to fit the low energy experimental data with 
the contribution from Coulomb force removed. Their nu- 
merical values are shown in table I. They are taken from 
Ref. (H: 



TABLE I: Scattering length a and effective range r for AW 
collisions. 



system 


parameter 


5 = 0,7 = 1 


S= 1,1 = 


PP 


a 


-17.1 (fm) 






r 


2.79 (fm) 




nn 


a 


-16.6 (fm) 






r 


2.84 (fm) 




np 


a 


-23.7 (fm) 


5.42 (fm) 




r 


2.73 (fm) 


1.73 (fm) 



In low energy scattering, the dominant contribution 
to the amplitude is given by the s-wave (orbital angular 
momentum £ = 0). Thus, we construct the spin-isospin 
averaged scattering amplitude from the expression (|24p 
for the s-wave with appropriate weight factors of spin and 
isospin. Details of the spin-isospin average are explained 
in Appendix A. 

In summary, the scattering amplitudes based on EFT 
are given by Eqs. (22]), ^23j), and (24]). It should be em- 
phasized that all these expressions of the scattering am- 
plitudes are valid only in a limited region of kinemat- 
ics. LO-ChPT (for tttt case) is usually considered to be 
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FIG. 1: (Color online) Comparison of the elastic cross sec- 
tions in 7T7T, ttN, and NN scatterings. Square and triangle 
points are the results of the low energy EFT and the phe- 
nomenological amplitudes, respectively. The actual fit to the 
experimental data was done for the differential cross sections. 
The NN elastic cross section is shown in logarithmic scale. 



Boltzmann equations will also be restricted to a small 
region in the T-fx plane, as we will discuss later. 

With this limitation in mind, it is quite important to 
consider physical scattering amplitudes in order to check 
the usefulness of, or more importantly to go beyond, 
the low energy EFT. To this end, we have constructed 
phenomenological amplitudes from the experimental data 
(differential cross sections) by fitting the coefficients of 
the partial-wave expansion. The fit was performed to the 
data up to y/s = 1.15 GeV (q cm — 550 MeV) in the 7T7T 
scattering, s/s = 2.00 GeV (q cm = 770 MeV) in the nN 
scattering, and y/s = 2.04 GeV (q cm = 405 MeV) in the 
NN scattering. Details of the fitting procedure are dis- 
cussed in Appendix A. In Fig.[TJ we compare elastic cross 
sections of the two different parametrizations. The elastic 
cross sections from the phenomenological amplitudes are 
almost identical with the experimental data (not shown) 
in the energy regions shown in the figures. Peaks in tttt 
and ttN cross sections are p-meson and A resonances, re- 
spectively. As mentioned above, the range of validity of 
LO-ChPT (for 7T7t, ttN cases) is p < 1 GeV for the pion 
momentum. Thus, if one applies this limitation to the 
momentum of colliding particles in the CM frame, one 
finds for the scattering energy \/s <C a/so with v /so ~ 2 
GeV for the 7T7t scattering, and y/so ~ 2.4 GeV for the 7rA^ 
scattering. Indeed, as evident from the figure, deviation 
of LO-ChPT from the phenomenological cross sections is 
already sizable well below the upper limits (mainly be- 
cause of the p-meson and A resonances). For example, 
in the 7T7t case, LO-ChPT gives the same tendency as 
the phenomenological cross section up to y/s ~ 600 MeV 
(9cm ~ 270 MeV), but beyond that the deviation is not 
small. Therefore, it would be safe to consider the validity 
region of LO-ChPT to be 



< 



s m FJ ee 600 MeV 



(7T7T scattering), (25) 



which is of course within y/s <C y/so ~ 2 GeV. Later we 
will use this limit to evaluate the maximum temperature 
up to which LO-ChPT is applicable. From the figures, 
it is obvious that the difference of two parametrizations 
becomes larger and larger with increasing energies. The 
7T7T and irN cross sections from LO-ChPT monotonically 
increase and become too large compared to the actual 
physical cross sections. We will see in the next section 
that this difference greatly affects the numerical value of 
the shear viscosity. 



valid for p <C 47r/„. ~ 1 GeV with p being the magni- 
tude of pion momentum. In the heavy baryon ChPT for 
the 7rA^ case, there is an additional expansion parameter 
p/rriN -C 1, which gives almost the same limitation as 
for the 7T7r case. Lastly, for the NN case, the differential 
cross section flMJ) is valid only for small q cm , or equiva- 
lently, small scattering energy near the threshold. This 
fact implies that the validity of results derived from the 



III. NUMERICAL RESULTS 
A. Pion gas 

Let us first discuss the case with only pions. This is 
important for properly understanding the effects of nu- 
cleons in the next subsection and, at the same time, for 
checking the validity of our calculation compared with 
the existing results [3,[iB|- As we mentioned before, our 
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formalism for the irN gas can be easily reduced to the 
case with only pions (by setting the nucleon degeneracy 
factor <?j\r - * 0, for example). Then, the shear viscosity 
coefficient is given by the first term of Eq. (fT5)) . and we 
determine B 7 ' (j>) by solving the Boltzmann equation ([l} 
with the collision term given by the first term of Eq. ([3]) . 
The expansion of B w (p) with respect to the orthogonal 
polynomials was taken up to the third order to ensure 
the convergence of the result. 



1. Range of validity 

Before presenting our numerical results, we clarify the 
range of validity in temperature for two parametrizations 
of the scattering amplitudes. As we already specified in 
Eq. (|25[) . LO-ChPT is valid only in a limited kinemat- 
ical regime: y/s ySmax = 600 MeV, while the phe- 
nomenological amplitude is by construction valid up to 
y/s = 1.15 GeV. These conditions can be translated to 
the limitation in temperature in the following way. Con- 
sider the binary collisions in thermal equilibrium. Since 
each collision takes place between particles in thermal 
equilibrium, the scattering energy squared s will fluctu- 
ate around its mean value (s) with a width S. The av- 
erage (s) and the width (standard deviation) £ may be 
defined by 



(s) 



(26) 
(27) 



Because the scattering energy squared of most of the 
collisions are below s max (T) = (s) + S, one can regard 
Smax(T) as (a measure of) the highest energy squared at 
temperature T. So, one may interpret the validity con- 
dition of LO-ChPT as s max (T) < s™£ = 0.36 GeV 2 , 
and thus obtain the validity condition for temperature 
T < T™? with given by s max (T^i^ ) = 0.36 GeV 2 . 

In Fig. [5J we show (s) and s max as functions of temper- 
ature. From the figure, one can read ~ 70 MeV. 
Therefore, we may conclude that the results of LO-ChPT 
are reliable only up to T ~ 70 MeV. On the other hand, 
the phenomenological amplitude is valid up to s ~ 1.3 
GeV 2 . Therefore, we expect from the figure that the re- 
sults of the phenomenological amplitude will be reliable 
up to temperature close to T c ~ 170 MeV. 

There is another important check for the validity of our 
framework. It is the applicability of the Boltzmann equa- 
tions. Recall that in deriving the Boltzmann equations, 
one assumes that a two-point correlation function can be 
factorized into a product of two one-point functions (one 
particle distributions f(x,p,t)). This is physically ac- 
ceptable when the density of particles n is small enough. 
This condition is normally expressed as 



<s> (GeV, 
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FIG. 2: (Color online) Temperature dependence of the aver- 
age scattering energy squared (s) and the measure of highest 
scattering energy squared s max = (s) + E. 
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FIG. 3: (Color online) The ratio A/rf as a function of temper- 
ature T. The description based on the Boltzmann equation 
is valid when A/c!> 1. 



where A is the mean-free path 



A = — 
no 



(29) 



A > d 



(28) 



with a being the cross section, and d is the interac- 
tion range, meaning that each collision takes place in- 
dependently. If we define the interaction range d by the 
Compton length of pions: d ~ l/m w , then the valid- 
ity condition for the Boltzmann equation reads X/d ~ 
m^/no 3> 1. Alternatively, if we define d through the 
cross sec tion a s a ~ nd 2 , we obtain another expression 
X/d ~ y / ir/(T 3 /n ^> 1 . Both expressions give similar re- 
striction on temperature. Figure [3] shows the ratio X/d 
as a function of temperature, where the (phenomenolog- 
ical) cross section a is estimated by its thermal average. 
One finds X/d <; 3 at T = 140 MeV, which manages to 
satisfy the inequalities. At this temperature, the mean- 
free path of pions is estimated as A ~ 4 fm, which is 
consistent with the literature [23J . As the temperature is 
decreased, the mean-free path becomes longer, and the 
Boltzmann description gets better and better. 
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Combining these two results, we may conclude that 
the phenomenological analysis of a pion gas based on 
the Boltzmann equation will be valid up to temperature 
T ~ 140 MeV which is much higher than the LO-ChPT 
limit T ~ 70 MeV, but lower than the critical temper- 
ature T c ~ 170 MeV. However, we notice that there is 
a caveat to this conclusion. In fact, even though we 
have satisfactory description of tttt scatterings in wider 
kinematical regime, we have to worry about at least two 
other effects as temperature increases. The first one is 
the effects of other (heavier) degrees of freedom such as 
kaons, and the second is the possible modification of pi- 
ons in thermal enviornmcnt. Both are however beyond 
the scope of the paper and we leave them for future prob- 
lems. 



2. Shear viscosity coefficient 77 

Figure 0] shows the shear viscosity coefficient 77 as a 
function of temperature T . Open diamonds (O) and open 
triangles (A) are based on numerical solutions to the 
Boltzmann equations with the phenomenological ampli- 
tudes and LO-ChPT, respectively. Remarkably, the two 
results show quite different behavior. While the result of 
LO-ChPT decreases with increasing T, that of the phe- 
nomenological amplitude shows the opposite behavior. 

Qualitative behavior of two different results can be 
easily understood by using rough estimate of the vis- 
cosity coefficient. In classical transport theory for a di- 
lute gas of one component, the shear viscosity coefficient 
is expressed in terms of the mean-free path A = 1/na, 
Eq. (HU), as follows: 



V 



-npX , 



(30) 



where n is the particle number density, p is the aver- 
age momentum, and a is the binary cross section. Thus, 
the shear viscosity coefficient is inversely proportional to 
the cross section. As we already saw in Fig. [2] the typ- 
ical energy in the tttt scattering increases with tempera- 
ture. This means that average cross section (a) indirectly 
depends on temperature. Notice also that the average 
momentum p = J d 3 p |p| fo(p)/ J (Ppfoip) roughly in- 
creases like p oc y/T (because (p) 2 /2m ~ 3/cT/2). There- 
fore, temperature dependence of the shear viscosity coef- 
ficient is essentially determined by the interplay between 
those of p and (a). For example, if the average cross 
section (a) increases rapidly as a function of T (faster 
than y/T), r\ ~ p] (a) is a decreasing function of T. But 
if (a) is almost constant, r\ is an increasing function of 
T. Based on these considerations and Fig. [TJ we can 
easily deduce the followings: As for the temperature de- 
pendence of the shear viscosity, we expect that ?7phcno is 
slightly smaller than ?7chPT at low temperature (because 
Cphono Z crchPT), while ?7phono > ?7chPT at high temper- 
ature. More precisely, since (a) of the ChPT monoton- 
ically increases as (a) ~ T 2 / f% while that of the phe- 
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FIG. 4: (Color online) Shear viscosity coefficient 77 of a pion 
function of temperature. Numerical results from the 
Boltzmann equations are compared with the rough estimate 
(ESI). 



nomenological amplitude does not grow, we expect that 
the shear viscosity decreases in the ChPT case while the 
opposite happens in the phenomenological case. 

Now let us come back to Fig. UJ where we also show the 
rough estimate (|30[) with two different parametrizations. 
For the (average) cross section a in the rough estimate 
(|30p . we evaluate it at the average momentum p, i.e., 
(a) = cr(p(r)).[32T] Comparing the results of the Boltz- 
mann equations and of the rough estimate (j30|) , we find 
that qualitative agreement of the results in two different 
ways of computation. Nontrivial behavior of the rough 
estimate for the phenomenological amplitude (filled dia- 
mond points) is due to the resonance shape of the elas- 
tic cross section in Fig. [TJ Indeed, the valley around 
T ~ 110 MeV corresponds to the p meson peak around 
s ~ 0.5 GeV 2 in the 7T7t cross section (According to Fig. [21 
s - 0.5 GeV 2 is translated into T - 150 MeV, which is 
further diminished due to p ~ y/T in the numerator of 
77). In the Boltzmann equations, such structure is further 
washed out by thermal average. From these analyses, we 
now understand that the behavior of the numerical re- 
sults of the Boltzmann equations is largely due to the 
energy dependence of the cross sections used in the colli- 
sion terms. In particular, the decreasing 77 of LO-ChPT 
at relatively high temperature T <; 80 MeV is an ar- 
tifact of too large cross section outside of the validity 
region of LO-ChPT. On the other hand, our most reli- 
able result (open diamonds O in Fig. [4j shows linear in- 
crease with temperature for T > 1 20 MeV. This behavior 
is consistent with that of hadronic resonance gas mod- 
els [1, Moreover, our result is also consistent with 
that of Ultrarelativistic Quantum Molecular Dynamics 
(UrQMD) for a meson gas [13], which also gives linear de- 
pendence on T, and rj ~ 0.1 GeV-fm" 2 = 0.0039 GeV" 3 
at T = 150 - 160 MeV. 



9 



EFT 

phenomenological 
KSS bound 



■o-o 



80 



100 



120 140 
T(MeV) 



160 



180 



FIG. 5: (Color online) The ratio rj/s as a function of tem- 
perature T in a pion gas system. Solid (A) and dashed (O) 
curves are, respectively, the results with LO-ChPT and the 
phenomenological amplitude. Dotted line corresponds to the 
conjectured lower bound rj/s = 1/471". 



3. The ratio rj/s 

Figure [5] shows our numerical results of the ratio 77/s as 
a function of temperature T. Dashed and solid lines cor- 
respond to the results of LO-ChPT and the phenomeno- 
logical amplitude, respectively. In both cases, rj/s is a 
monotonically decreasing function of T. Notice that the 
result of LO-ChPT violates the conjectured bound 1/4-77 
(the KSS bound, shown as the dotted line) at around 
the critical temperature T c ~ 170 MeV. On the other 
hand, 77/s from the phenomenological amplitude keeps 
well above the KSS bound up to temperature ~ T c . 
These are consistent with the results of Refs. (tS. Hij. 

Since the entropy density is common in both cases s oc 
T 3 , behavior of 77/s can be understood by that of 77 itself. 
For example, we saw that 77 decreases in LO-ChPT, while 
it increases in the phenomenological case. Such difference 
affects on the ratio rj/s: it decreases faster in the LO- 
ChPT case than in the phenomenological one. Also there 
is a crossing point for the two results of 77/s at T ~ 
80 MeV (see Fig. [5]), and this point coincides with that 
of the shear viscosity 77 (see the open diamonds and open 
triangles in Fig. [3J . 



It has been argued in Ref. [14J that violation of the KSS 
bound 77/s > 1/47T suggests the existence of phase transi- 
tion (or crossover transition) in order for the KSS bound 
to remain valid. This kind of argument is of course dan- 
gerous because the precise value of 77/s depends on the 
amplitude in the collision term, and the result of Ref. [lij 
is based on LO-ChPT, whose applicability is limited to 
T ^ 70 MeV as we already discussed in detail. In fact, 
more reliable result with the phenomenological amplitude 
does not violate the KSS bound even around the critical 
temperature. Therefore, "violation of the KSS bound" 
cannot be the signature of phase transition. Instead of 
seeing violation of the KSS bound, we propose to check 
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FIG. 6: (Color online) The absolute value of the slope 
\d(jj I s) I dT\ for the result of phenomenological amplitude. 
Straight line is the linear extrapolation at T = 140 MeV. 



the other property to catch the indication of phase tran- 
sition. As we mentioned in Introduction, Ref. [6( suggests 
that the ratio has a minimum around the (pseudo-) crit- 
ical temperature T c . If this is indeed the case, it implies 
that (for crossover transition) the curve of 77/s as a func- 
tion of temperature will have a convex form around T c , 
and the slope will decrease as temperature approaches T c 
from the left. (This kind of argument will hold as far as 
the temperature is not too far away from the critical end 
point.) What we observed in Fig. [5] for the phenomeno- 
logical amplitude is indeed the decrease of the slope with 
increasing temperature. Thus, even if the ratio is still 
well above the conjectured bound, we can anticipate the 
existence of phase transition. 

In Fig. [6l we show the (absolute value of) slope of 
the ratio 77/s for the phenomenological amplitude. Re- 
call that we have estimated the border of validity region 
of our calculation to be T ~ 140 MeV. Thus, we use 
the results around T ~ 140 MeV for linear extrapolation 
towards higher temperature. More precisely, we approxi- 
mate the curve 1Z(T) = rj/s around at some temperature 
T = T < T c as 



U(T) ~ n(T )+n.'(T )(T-T ) + 



11" (T ) 



(T-T ) 2 . (31) 



Then the slope dlZ(T)/dT is approximated by a lin- 
ear function of T. The critical temperature T c may 
be defined by the temperature where the slope is zero: 
dK(T)/dT = 0. Namely, 



T r . ~ Tn 



K'(Tp) 

n"{T y 



(32) 



where W(T ) < and TZ"(T Q ) > on the left-hand side 
of a convex function 1Z(T). The result of linear extrapo- 
lation at T = To = 140 MeV is shown on the same figure. 
Remarkably, the temperature at which the straight line 
cuts the horizontal axis is 173 MeV, which is quite a rea- 
sonable result as the critical temperature. In addition to 
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this, one can guess the value of rj/s at T = T c by using 
the approximation (f3Tj) : 

K(T C )^K(T )- ^^ . (33) 

If one substitutes numerical values at To = 140 MeV, 
namely, K(T ) ~ 1.0, W(T ) ~ -0.005 MeV" 1 and 
ft"(T ) ~ 1.38 x 10~ 4 MeV" 2 which are read from Figs.E] 
and[5J one finds K = rj/s ~ 0.9 at T = T c . These re- 
sults should be understood with reservation at least for 
two reasons. First of all, there is ambiguity in the choice 
of To, which will affect the values of T c and rj/s. How- 
ever, in fact, T) cannot be taken arbitrary because we 
need to take To as close to T c as possible for the linear 
extrapolation to be accurate. To obtain reasonable val- 
ues of T c and rj/s, the largest possible value for T is 
preferable. Since the choice To =140 MeV is the upper 
limit of our validity region, we expect that the estimated 
value T c = 173 MeV is the best value of our calcula- 
tion. Therefore, even if there might be some ambiguity 
in selecting To, we can say that we have chosen the best 
value. The second source which might change the values 
of T c and rj/s is the possible contributions from heavier 
mesons. Since such contributions become more impor- 
tant as T — > T c , our results with only pions become bet- 
ter as we depart from T c (which is however not desirable 
for determination of T c ). Still, we expect our results are 
not so bad because at T ~ 100 MeV such heavier parti- 
cles can be ignored, and even at T = 140 MeV, number 
of kaons amounts to only 20% of total particles. 

For the crossover transition, the ratio rj/s will be con- 
tinuous at the (pseudo) critical temperature T c . This 
immediately implies that the numerical value of rj/s de- 
termined above is relevant even in the deconfined phase. 
Our result rj/s ~ 0.9 is well above the KSS bound, but 
is small enough compared to the weak-coupling QCD re- 
sult (see for example, Fig. 4 of Ref. jfj}). In this sense, 
the QCD matter around T c could be "strongly interact- 
ing" . However, we should be careful when we draw such 
a conclusion from the value of rj/s. In fact, the small- 
ness of the ratio n/s is not a direct consequence of large 
cross section which may be realized by a strongly inter- 
acting matter. In our calculation with the phenomeno- 
logical amplitude, the cross section does not grow a lot 
(in contrast with the LO-ChPT) and the viscosity r\ even 
increases as T — * T c . Still, since the entropy increases 
faster than 77, the ratio rj/s becomes a decreasing func- 
tion of temperature. Therefore, the smallness of the ratio 
is realized in a nontrivial way. 



B. Pion-nucleon gas 

Let us now present the numerical results for a dilute 
7riV gas. As advocated in Introduction, addition of nu- 
cleons to a pion gas enables us to study the effects of 
baryon chemical potential. We discuss how the results of 
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FIG. 7: A measure of the highest scattering energy squared 
= (s NN ) + Z NN for the NN scattering is plotted in 
the T-jj, plane. Gradation represents the values in the range 
3.5 GeV 2 < <7 GeV 2 . Solid curve corresponds to the 

borderline s££ = 5.29 GeV 2 . 



a pion gas presented in the previous subsection change 
under the influence of the chemical potential. 

1. Range of validity 

We start again by the discussion on the range of valid- 
ity of our framework. Since we have to take the ttN and 
NN collisions into account in the 7riV gas, we introduce 
the average scattering energy squared between particles 
i and j — ir or N): 

s J^J^^(Pi,ft)fo(Pi)/^) 

For the -k-k scattering, this is of course equivalent to 
Eq. (|26[) and depends only on T. But for the ttN and 
NN scatterings, the average values depend on both T 
and ji. Similarly as before, we further define the stan- 
dard deviation by E y = yj ((s y ) 2 ) — (s' 1J ) 2 , and regard 
s^ ax (T, (i) = (s u ) + E y as a measure of the highest en- 
ergy of the i,j scattering at temperature T and baryon 
chemical potential /1. 

In Fig. [3 we show the values of s^. on the T-ji 
plane. Recall that the fit to experimental data in the 
NN scatterings is by construction valid up to s NN < 
(2.04) 2 = 4.16 GeV 2 . However, in fact, our parametriza- 
tion works well up to slightly higher value s NN ~ (2.3) 2 = 
5.29 GeV 2 . Therefore, we define the borderline of the va- 
lidity region by temperature and chemical potential that 
satisfy s££(T,/i) = 5.29 GeV 2 . The result is shown on 
the same figure [7] by a thick solid curve. The maximum 
temperature is about 130 MeV or slightly smaller than 
that, while the chemical potential is not restricted up to 
l-i = l GeV. If wc perform the same analysis for the ttN 
scattering, the borderline defined by sj^ x = 4 GeV 2 lo- 
cates outside of the region of our interest T < 180 MeV, 
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FIG. 8: (Color online) The lines for the ratio A^jv/d in a irN 
gas. In the region of low density and high temperature, we 
show only X n /d, while in the region of low temperature and 
high density, Ajv/d. We may regard the line for X/d — 3 as 
the border of the validity region of the Boltzmann equations. 
The line for X^/d = 1 is far outside of the T-y, region shown 
here. 
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FIG. 9: (Color online) Temperature dependence of the shear 
viscosity coefficient r\ of a irN gas at different values of chemi- 
cal potential /i = 300, 500, 700 MeV. Comparison is made be- 
tween the results of the low energy EFT and the phenomeno- 
logical amplitudes. 



jit < 1 GcV. On the other hand, the validity region of 
the low energy EFT is severely restricted. For example, 
if we take the maximum value of V s NN to be 1.88 GeV 
(beyond which the phenomcnological fit starts to devi- 
ate from the low energy EFT fit, see Appendix A), then 
the borderline defined by s££ = (1.88) 2 = 3.53 GeV 2 
is very close to the horizontal axis (see Fig. [7]) . Even if 
we relax the condition to higher value = (1.90) 2 = 
3.61 GeV 2 , allowed region is still very narrow. 

Let us also examine the validity condition for the Boltz- 
mann equations. The criterion is again given by Eq. (|28[) . 
but now pion's mean-free path should be modified in the 
presence of nucleons, and we have to separately consider 
nucleon's mean-free path, too. According to classical 
transport theories, the mean-free path of i-th component 
in a mixed gas is modified as [25[ 



A, 



A? 



1 + £i* 



aij rij 
2 Uii m 



(35) 



where A? = l/nj<Tj is the mean- free path of a pure gas of 
i-th component, and is the cross section between i-th 
and j-th components. Therefore, the mean-free paths of 
pions and nucleons are respectively given by 



1 



, l+m^ /tun 

n n a n7T + ting^n \l f 



VV 



1+mjv /m n 

nN<?NN + n v a vN \/ ^ 



(36) 
(37) 



Consider the condition for nucleons. If one takes 
cat at ~ 40 mb as a typical value of the NN cross section 



(twice of the saturating value at high energy, see Fig. [T]) 
and uses d ~ l/m ff again for the interaction range, then 
the condition for nucleons A at 3> d can be estimated as 
n 7v <C rn-n/oNN — n% with = 0.157 fm -3 being the 
normal nuclear density. This condition is physically quite 
reasonable since we do not expect standard Boltzmann 
description useful at normal nuclear matter density. To 
obtain more precise restriction depending on T and /i, 
we need to use Eqs. (f3"6"f and (|37[) and estimate each 
condition by replacing a and n by their thermal aver- 
ages. In Fig. we show the lines for several values of 
the ratio X/d in the T-fi plane. For simplicity, we used 
d = 1/™^ for both pions and nucleons. In the region of 
low baryon density where pions are dominant degrees of 
freedom for transport phenomena, we show the ratio for 
pions X^/d. The line for X^/d = 3 (solid line) is consis- 
tent with the previous result shown in Fig. [3J We may 
regard this line as the border of the validity of the Boltz- 
mann equations at low densities. On the other hand, at 
large chemical potential and low temperature, transport 
phenomena is dominated by nucleons. The lines for the 
ratio Xff/d — 1, 3, 10 are shown in the figure. We may 
again regard the line Xpf/d= 3 as the borderline, which 
reaches at (T, /x) = (0, 950MeV). This condition is more 
restrictive than that of Fig. [7] 

Combining these two analyses, we may conclude that 
our framework is valid in a wide region of the T-fi plane, 
whose boundary is roughly given by (a quarter of) the 
elliptic curve connecting (T, fi) ~ (130MeV, 0) and (0, 
950MeV). The first point (130MeV, 0) is specified by the 
measure of the highest scattering energy squared, and the 
last point (0, 950MeV) is from the validity limit of the 
Boltzmann equations. 
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FIG. 10: (Color online) fi dependence ol 77 at T = 50MeV (left) and T = lOOMeV (right), and its decomposition 77 = -rf + 77' 



2. Shear viscosity coefficient 77 

In Fig. [51 we show temperature dependence of 77 at 
different values of chemical potential li = 300, 500, 700 
MeV. Comparison is made between the results of the 
low energy EFT and the phenomenological amplitudes. 
These temperature dependence is qualitatively consistent 
with the previous results of the pion gas (see Fig. [4|. As 
for the /i dependence, however, the shear viscosity coef- 
ficient shows nontrivial behavior. In particular, in the 
window 80 MeV < T <130 MeV, it increases with in- 
creasing ll. This is not quite understandable at first be- 
cause we expect that the cross section of tttt scattering 
will effectively enhance in the presence of nucleons while 
the effects of nucleon viscosity may be ignored at lower 
density, meaning that the shear viscosity will decrease. 
To understand what really happens when \i =^ 0, let us 
again consider a rough estimate of the shear viscosity of 
a mixed gas. Let the shear viscosity coefficient of a pure 
gas of particle species % be t/ = riip^/3 where n, is the 
number density, p i is the average momentum, and A° is 
the mean-free path. Then the shear viscosity coefficient 
for n component gas is given by the sum of each visc osity 
t]q with modified mean- free path given in Eq. (|35p [25j : 

i 1 

Thus for the ttN gas mixture, we obtain 

r) = rf + r) N (38) 
Vo , rtf ^ 

1 j 1 ( OxW ) ( n N ) 1 1 / ni N I A I n„ \ 

where we have used the approximation m^/mN *C 1. In 
the two extreme limits n^jn^ — > and 00, the formula 
(|38p reduces to 77J and 77^, respectively. Therefore, this 



formula interpolates a pure pion gas at low fi and a pure 
nucleon gas at high /i. Notice that the pion contribution 
rf is always smaller than 77J due to the presence of nN 
interaction, as we alluded before in relation to Eq. (|2"Tj) . 
This is exactly what we expected. On the other hand, 
it is not straightforward to predict the behavior of the 
total shear viscosity. If one plots the rough estimate (f38|) 
as a function of un/ti^ assuming that the cross sections 
are constant and are of the same order, one finds that 77 
decreases at small values of jn^ (small /z), but turns 
into increase at large un/h^ (large fi). If one changes the 
numerical value of cross sections a little, then the curve 
turns into monotonic increase. In fact, both can happen 
in reality depending on temperature. In Fig. 1 1 01 we have 
plotted ry at T = 50 MeV and 100 MeV as a function of 
\x. While the total viscosity behaves differently in these 
two panels, one can see monotonic decrease (increase) 
of rf (r) N ), namely, the interplay between rj 11 and 77^ 
at both temperatures. Therefore, the increase of 77 with 
increasing [i observed in Fig. [5] can be understood as a 
result of enhancement of 77^ . 



3. The ratio rj/s 

Figure [11] shows temperature dependence of the ratio 
77/s of the irN gas at different values of baryon chemi- 
cal potential ll = 300, 500, 700 MeV. The left-hand side 
is the results of the low energy EFT, while the right- 
hand side, the phenomenological amplitudes. The ratio 
77/s shows qualitatively the same behavior as in the pion 
gas system (Fig. [5]): (i) Both the results (EFT and phe- 
nomenological amplitudes) are monotonically decreasing 
functions of T, (ii) (?7/s)eft > (»7/s)phcno at lower T 
while opposite at higher T, and (iii) the ratio of the 
EFT violates the KSS bound at around T ~ T c while 
that of the phenomenological amplitudes does not. No- 
tice that the inclusion of chemical potential works to re- 
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FIG. 11: (Color online) The ratio n/s of a nN system plotted as a function of temperature T at different values of baryon 
chemical potential li = 300, 500, 700 MeV. (Left) results of low energy EFT, (Right) results of phenomenological amplitudes. 



duce the value of rj/s. As a result, for the ratio of the 
EFT, the temperature at which the curve cuts the KSS 
bound becomes smaller. On the other hand, for the re- 
sult of the phenomenological amplitudes, the flattening 
of the curves seems to occur at lower temperature with 
increasing /i. If this is indeed the symptom of crossover 
transition as we discussed before, one can say that the 
(pseudo) critical temperature T c will decrease with in- 
creasing fx, which is consistent with what we know from 
lattice simulations or effective models. [33j 

Our result (?7/s)p hono ~ 0.5 - 0.4 at li = 700MeV 
and T ~ 100 - 140 MeV is consistent with that of 
"URASiMA" (Ultra-Relativistic AA collision Simulator 
based on Multiple scattering Algorithm) [Lsj which is a 
Monte-Carlo event generator of hadronic collisions. It in- 
cludes both elastic and inelastic scatterings whose cross 
sections are given by experimental data. The numerical 
coincidence of the ratio from different frameworks is very 
interesting. In fact, the essential difference of URASiMA 
from our framework is the presence of inelastic collisions. 
But as far as we consider small deviation from thermal 
equilibrium, inelastic collisions which will change particle 
numbers (such as rnr ~ * 7T7r7r or ttN — ► ttttN) would not 
be so important, and it seems reasonable to obtain the 
same result from two different frameworks. 

As we already saw in Fig. [91 the shear viscosity coef- 
ficient increases in the window 80 MeV< T <130 MeV 
with increasing li. However, in the same window, the ra- 
tio does decrease. This clearly implies that the reduction 
of rj/s at T ~ 80 — 130 MeV is due to the increase of 
entropy. 

In Fig. [T21 we show the /i dependence of the ratio rj/s 
at temperature T — 50, 100 MeV (left) and at low tem- 
perature T = 10 MeV (right). Recall that the ratio de- 
creases with increasing li in Fig. QT] where T > 70 MeV 
was shown. This is consistent with the left figure, and 
in agreement with our original expectation as we men- 
tioned in Introduction. On the other hand, the right 



figure shows a new structure: There is a valley at large 
Li ~ 950 MeV. We will discuss later the physical implica- 
tion of this structure. 

Let us comment again the point made in the last para- 
graph of the previous subsection. The ratio rj/s becomes 
less than 0.3 at T > 140 MeV and fi = 700 MeV. 
This is small enough and is close to the KSS bound 
rj/s ~ 0.1 compared to other systems such as water. 
However, as shown in Figs. [5] and [TU1 the shear viscos- 
ity itself grows as the system approaches phase boundary 
(T — > T c , li — > Li c )- Therefore, even if the ratio is small 
enough, the 7rV system cannot be treated, strictly speak- 
ing, as a perfect fluid with r\ — 0. Sometimes ideal hydro- 
dynamics is used to describe the matter after hadroniza- 
tion in heavy- ion collisions, but one will have to take into 
account the effects of viscosity for a realistic simulation. 



4- Valley structure at large /i and low T 

In the right panel of Fig. [T^] (low temperature 
T = 10 MeV), we pointed out a valley structure at 
high baryon chemical potential, which was not seen in 
higher temperature (left panel). Let us look at this 
new structure in more detail and examine its possible 
interpretation. We recall again Ref. @, where it was 
suggested that the ratio rj/s would give a minimum 
at the phase transition temperature. If this is true 
for other control parameters, in particular, chemical 
potential, and if there is phase transition under the 
change of chemical potential, it is natural to expect 
that the ratio would exhibit a valley structure with 
its minimum at the critical chemical potential ji c . In 
other words, if one finds a valley structure in the yi 
direction, one can expect some kind of phase transition 
around the minimum. This is what we observed in 
Fig. [T2] Then, what kind of phase transition could be 
related to this valley structure? The valley locates at 
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FIG. 12: (Color online) [s, dependence of n/s: (left): the results of the phenomenological amplitudes at higher temperature 
T = 50, 100 MeV, (right): the results at lower temperature T = 10 MeV. 



low temperature T < 20 MeV and at high chemical 
potential /1 ~ 950 MeV (which is however not enough for 
the quark- hadron phase transition). This is the re gion 
where we can see the liquid-gas phase transition. [34| 
At low T < 20 MeV and around normal nuclear density, 
there is a critical line (possibly first order) separating a 
nucleon gas phase and a nuclear matter (liquid) . This 
line terminates at around T ~ 15 MeV [23], and above 
that temperature, there is no distinction between a gas 
and a liquid. Therefore, if the valley indeed corresponds 
to the liquid-gas phase transition, it should disappear 
when the temperature go far beyond T ~ 15 MeV. In 
Fig. [13l we show the transition of the valley structure 
from T = 5 MeV up to T = 20 MeV. The right panel 
is the results of the phenomenological amplitudes. 
Clearly, with increasing temperature, the valley becomes 
shallow, which supports the interpretation that the 
valley structure indeed corresponds to the liquid-gas 
phase transition. There are several comments about this. 

• Although the liquid-gas phase transition is a phe- 
nomenon in the hadronic phase, it is not obvious 
whether the Boltzmann equations (valid for a dilute 
gas) correctly describe the transition to liquid phase. 
Notice that the region where liquid-gas phase transition 
takes place is close to the border of the validity region of 
the Boltzmann equations, and thus it is not surprising 
that there might exist another phase (i.e., liquid phase) 
outside the region of validity. Still, it would be safe 
to reserve that reliable information from our analysis 
should be only the tendency towards phase transition. 
We do not expect we can describe precise structure of 
phase transition such as the order of transition and 
the temperature dependence of the critical chemical 
potential /x c (T) (Our result is that the minimum of the 
valley moves to the right with increasing T, as opposed 
to common expectation). In order to describe the 



transition correctly, we will have to include the effects 
of higher-order correlations to the Boltzmann equations, 
or start from different models (such as the a-cu model) 
which are more appropriate for nuclear matter, both 
of which, however, are beyond the scope of the present 
paper. 

• As mentioned in Introduction, a similar valley 
structure was already reported in the calculation with 
the low energy EFT [l6j | and the authors of Ref. [lj| 
claim that it is related to the liquid-gas phase tran- 
sition. We have done the same calculation with our 
own parametrizations of low energy EFT, and obtained 
consistent results with Ref. [lj| as shown in the left panel 
of Fig. [13J Unlike the results of the phenomenological 
amplitudes (the right panel), the valley persists even 
at higher temperature. As we repeated several times, 
since the range of validity of the low energy EFT is 
severely restricted (see the discussion about Fig. [7]), it 
is quite dangerous to draw any conclusions about the 
physics outside of the validity region. However, a similar 
valley structure is observed even in our calculation with 
the phenomenological amplitudes, and thus our cal- 
culation partially supports the results of low energy EFT. 

• Since the entropy in our calculation is evaluated by 
using free particle distributions fn' (p), there is no in- 
formation about phase transition in the denominator of 
r)/s. However, if there is a real phase transition, entropy 
will of course show nontrivial change around T c or [i c . 
For example, in the liquid-gas transition of water, both 
the shear viscosity and the entropy have nontrivial struc- 
ture around T c , and contribute to give a convex shape of 
the ratio 77/ ' s. We emphasize that, in our calculation, the 
valley structure of n/s is not the result of entropy. If 
we treated actual entropy, the structure would emerge in 
more pronounced way. 




FIG. 13: (Color online) Temperature dependence of the valley structure of rj/s at relatively high chemical potential. Only the 
results at low temperature T = 5, 10, 15, 20 Me V are shown. (Left): Low energy EFT, (Right): Phenomenological amplitudes. 



IV. SUMMARY 

We have performed a detailed calculation of the shear 
viscosity coefficient r\ and the viscosity to entropy ratio 
77/s in a wide region of the hadronic T-fi plane. Our for- 
malism is based on the relativistic quantum Boltzmann 
equations, and we found that it is very important to use 
phenomenological amplitudes for the scattering ampli- 
tudes in the collision terms in order to obtain reliable 
results. On the other hand, the validity region of the low 
energy effective field theories is severely restricted, and 
reliable results based on them are also limited to a small 
region of T-fi plane. We found that the ratio r\j s de- 
creases (for T ?Z 20 MeV) under the inclusion of nucleon 
degrees of freedom, but still keeps above the conjectured 
KSS bound r]/s — l/4ir in the region we investigated 
(T < 180 MeV, /i < 1000 MeV). Since the shear viscosity 
coefficient itself increases with increasing temperature, 
the behavior of the ratio is largely due to the entropy. 
At low temperature T < 15 MeV and high baryon chem- 
ical potential fi ~ 950 MeV, we found a valley structure 
in the ratio. There is some argument about the relation- 
ship between such structure and phase transition, and we 
expect that the valley structure found in our calculation 
would correspond to the liquid-gas phase transition. 
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APPENDIX A: PHENOMENOLOGICAL 
AMPLITUDES 

In this Appendix, we explain the phenomenological 
amplitudes we used in the kinetic equations. Since the 
scattering amplitude (squared) is related to the differen- 
tial cross section in the CM frame as 

with s being the Mandelstam variable (scattering energy 
squared), we discuss only the differential cross sections. 



1. 7T7T scattering 
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Consider the partial- wave expansion of the isospin av- 
eraged 7T7T elastic differential cross section. Taking up to 
p-wave scattering (£ — 1) yields a very nice description 
of the experimental data: 
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(19. 



2 

==— — — £(2/+l) — £ |(2^ + l)e^sin^P,(cos^ cm ) 

li>k £1 + 1 ; J=0 1™ 1=0,1 

sin 2 6° Q + J sin 2 «5 2 + | • 9 sin 2 # cos 2 cm ) , (A2) 

9cm \ 9 9 J / 



where 4 in the overall factor comes from the identical fac- 
tor, cm and q cm are the scattering angle and the magni- 
tude of momentum in the CM frame, and Pi is the £-th 
order Legendre polynomial (Po(x) = 1, P\{x) = x). The 
numerical factors in the last line are from isospin / and 
orbital angular momentum i. For example, the last term 
corresponds to J = I = 1 scattering and thus (3/9) • 9 
is from the isospin weight (21 + 1)/ Y<r( 21 ' + 1) = 3/9 
for / = 1 and (2£ + l) 2 = 9 for I = 1. The phase shift 
8\ depends on i and /. Since the total wavefunction of a 
7T7T system must be symmetric under the exchange, only 
three phase shifts (I, t) = (0, 0), (1, 1), (2, 0) are possible. 
It is known that the energy dependence of the phase shifts 
are well parametrized by the following function [28j : 



ential cross section: 



tan 5l 



- \l 1 — 9cm 



Ami 



(A\ + B\q 



2 

cm 



t 9cm 



s - si 



(A3) 



where the parameters are determined to fit the data up 
to y/s = 1.15 GeV (<7 cm ~ 550 MeV) and are shown in 
table II: 



TABLE II: 





{I, l) = (0,0) 




(J,*) = (1,1) 




(!,£) = (2,0) 


A 


2.25 x KT 1 


Ai 


3.63 x lGT* 


Ai 


-3.71 x 10~ 2 


B° 


2.46 x 10 _1 


B{ 


1.34 x 10~ 4 


Bl 


-8.55 x 10 -2 


C° 


-1.67 x 1(T 2 


cl 


-6.98 x 10~ 5 


cl 


-7.54 x 10" 3 


D° 


-6.40 x 1(T 4 


Dl 


1.41 x 10~ 6 


Dl 


-1.99 x 10" 4 


si! 


36.7 


si 


30.7 


si 


-11.9 



In this table, dimensionful parameters are redefined so 
that the mass dimension is provided by m n . For ex- 
ample, since Bq has mass dimension —2, we define 



b; 



o = b° /r 



i 2 and 6q = 2.46 x 10 _1 . This phenomeno- 
logical parametrization describes the experimental data 
very well. We use this for the irir scattering amplitude. 



2. iv N scattering 

Let us now turn to the irN scattering. Consider the 
partial- wave expansion of the isospin averaged irN differ- 



da_y 

dSl J 



N 



Hem £=Q 



-.1=1/2 



(9c 



xP^(cos( 



The coefficients 



Ci{q c 



= gQ£ /2 (9c 



+ gQ' /2 (9c 



ft (9cm) 

(A4) 





are functions of q cm and are determined from the ex- 
perimental data. More precisely, each function Ci is ex- 
pressed by a superposition of 30 Gaussians: For example, 
for I = 0, we use 
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C!«=o(?cm) = ^ A n cxp I - 



71=1 

30 



gcm-lOOMeV V 
60"MeV ) 



n =w I V 60«MeV J /' 

where positions and widths of Gaussians are found by 
trial and error. We use different values of positions and 
widths for different £. This kind of technique is some- 
times used in describing nuclear many body wavefunc- 
tions. By using these functions and taking the maximum 
angular momentum 2£ maK = 8, we can fit the experimen- 
tal data up to Vs = 2 GeV {q cm ~ 770 MeV). Notice 
that this parametrization works very well at relatively 
high scattering energies, but in fact its quality becomes 
worse at small scattering energies. This is due to the fac- 
tor l/g 2 m in Eq. (|A4[) . If one absorbed this factor into 
the coefficients and performed the same Gaussian fitting, 
quality of the fit would be better even at low scattering 
energies. However, we decided to start from the conven- 
tional expression shown in Eq. (|A4|) , and to find another 
parametrization for the low energy data. We have inter- 
polated the parametrization proposed in Ref. [30j which 
is compactly represented for the scattering amplitude: 

M nN = b + h(t- f) + (c + ci (t • f) W • q') 
+i (do + di (t- T}ja-(qxq'), 

where t and f/2 are the isospin vectors of a pion and a 
nucleon respectively, q and q' are the momenta of incom- 
ing particles in the CM frame, and lastly a is the Pauli 
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TABLE III: 



bo 


-0.010 


bi 


-0.091 


CO 


0.208 


Cl 


0.175 


do 


-0.190 


di 


-0.069 



matrix. The parameters are determined as shown in ta- 
ble III (dimensionful parameters are again redefined by 
using so that they become dimensionless) . 
These two different parametrizations are smoothly 
matched at y/s = 1101 MeV (q cm =79 MeV), giving a 
very nice parametrization of the experimental data for a 
wide region of scattering energies. 



3. NN scattering 

Let us finally discuss the NN scattering. In the text 
we discussed only the s-wave scattering, but here we treat 
arbitrary orbital angular momentum. We define the spin- 
isospin averaged NN differential cross section 



da 



NN 



averagee 



~ 16 | [dQ + \dfl 



1,0 



o (da 

+3 U 



da 



,(A5) 



where 16 = £/=o,i( 2/ + l)£ ff=0 ,i(25 + x ) and 
[da / dVL) 1 ' S in the right-hand side are the differential 
cross sections with isospin / and spin S specified. One 
can further decompose each contribution depending on 
the value of orbital angular momentum £. (Notice that 
the total NN system must be antisymmetric under the 
exchange of two particles: (— \y+ s + I = — 1.) 
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dn y 
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Hera. 
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1,0 
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dfi, 
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0,1 
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dn y 


) - 




d<T S 


1,1 
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) - 


9cm 



£— even 



Inserting these into Eq. (IA5|) . we obtain the spin-isospin 
averaged differential cross section separately for even and 
odd I: 



do 

dn 



NN 



averaged 



r 



\{cos9) 



(A10) 



When I = 0, this gives the spin-isospin averaged scatter- 
ing amplitude used in the parametrization of low energy 
effective theory. 

Restriction to the s-wave alone may be a good ap- 
proximation at low scattering energy, but to obtain a 
parametrization which describes the data in much wider 
range of energies, we have to include larger I. To do this, 
we again perform the partial- wave expansion of the spin- 
isospin averaged differential cross section (IA5[) and fit the 
energy-dependent coefficients by using the Gaussian su- 
perposition. We define the coefficients Di(q CTn ) similarly 
to Eq. (fXH) : 

d \ NN 1 — « 

■£) = —J2 D dQcnr)Pi(cOs9). (All) 

aS2 / averaged 9cm g 

Since the low energy data are already described well by 
the parametrization (|24[) and the parameters given in ta- 
ble I, we adopt them as the low energy part of the global 
parametrization. We use them up to q cm — 6.76 MeV 
(y/s = 1.88 GeV), and beyond that, we switch to the 



partial- wave expansion (jAlip . For actual fitting of the 
experimental data [3l|, we divide the rest of the region 
into two: 

(i) 6.762 MeV < q cm < 48.76 MeV 
(1876 MeV < < 1879 MeV) 

(ii) 48.76 MeV < g cm < 405 MeV 
(1879 MeV < y/s < 2043 MeV) 

In each region, we determine the coefficient Di{q cm ) sim- 
ilarly as in the case for irN scattering, so that the curve is 
smoothly connected at the matching points q cra = 6.762 
MeV and <j cm = 48.76 MeV. In region (i), the highest 
value of the angular momentum is taken to be 6, while 
in region (ii), much higher value 16. We performed the 
fit up to q cm = 405 MeV (y/s = 2043 MeV). 



APPENDIX B: SOLVING THE LINEARIZED 
BOLTZMANN EQUATIONS 

In this Appendix, we discuss how to solve the Boltz- 
mann equations fT]) and @ . First of all, recall that the lo- 
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cal equilibrium state is defined by the distributions f£' 
which make the collision terms (the right hand side of 
the Boltzmann equations) vanishing [see Eqs. ©, (|7J)]. 
However it should be noticed that the local equilibrium 
is not the solution to the Boltzmann equations: Indeed, 
the parameters T, /i and characterizing f$' N are de- 
pendent on the coordinates, which implies that the left 



hand side of the Boltzmann equations (which have coor- 
dinate derivative d^) do not vanish. With this in mind, 
the Boltzmann equations linearized with respect to the 
deviations 8f v ' N from the local equilibrium f^' N (x,p) 
can read as follows (leading order of the left hand side 
does not contain 6f*' N ): 



= c™[sr, /ol + c^LfoV/l + c* N [sr, f N ] + c N [fS, *f N ] > < B 1 ■> 



j^dX = C N *[6f N , f*] + C N *[f N , &n + C NN [Sf N , f N ] + C NN [f N , Sf N ] , (B2) 
v 

where the notation for the collision term is for example (see Eq. ||3J)) 

c™[sr,f5} = c™ifz+6rjz]-c™[fsjz) 

+ fa)f&f& - (i + /oi)(l + fSt)SfSfSp} ■ (B3) 

This corresponds to the lowest order Chapman-Enskog method. So far, the deviations 5f' are in principle arbitrary, 
but for the purpose of computing the shear viscosity coefficient rj, we can restrict only to the deviations that are directly 
from the shear diVj ^ (i ^ j). Thus we ignore any other effects except the shear. 

By using the particle number conservation dtn = — iiViV % and the energy-momentum conservation dJT^" = 0, the 
left hand sides of Eqs. (|B1[) and (|B2|) can be rewritten as 



f^«(x) = f -^±^ L iPj - | P 2 ) (V«^) trl , ( B4 ) 



jgOXM - P f ° N{1 EpN f ° N) (w - 6 f P 2 ) ( W)« > ( B5 ) 

where (V'V 3 '). , is defined in Eq. pill . Notice that this functional form was the motivation for defining new quantities 
B 7T N {p) in Eqs. (fl3 ]> . (jig]) . Indeed, with respect to B n > N (p), the right hand sides of Eqs. (jBTJ and (|B2) can be expressed 
rather compactly. Introducing further the following notation, 

,N, \ _ ( - - &ij 



B% N (p)^(p i p j - u f)B«> N (p), (B6) 
the linearized Boltzmann equations (|B 1 1) and {B2| can be expressed as (/^ = fo(k\), f$ p = fo(p), etc.) 

fSp{1 ^ fSp) (p iPj - S fp^ = i^J dT*(l + + f£ 2 )fS 3 f£ p (B*(p) + B-(k 3 ) - £g.(fe) - B%(kij) 

+ g N f dT« N {l - /£)(! + fS 2 )f? 3 fS p (B%(p) + B% (k 3 ) - B*(k 2 ) - B% (k x )) ,(B7) 



f ° Npil E J° p) (piPj- 5 fp 2 ) = 9* J ^*(l + /oW-/S)/S/5( B «(P) + ^(*3)-^(fe)- fl 5(*i) 

v (l - /£)(! - (i$ (p) + B? (fe) - S£ (fc 2 ) - B£ (fex)) .(B8) 



H / dT NN < 



These are the equations for B n ' N (p). As explained in to a finite dimensional functional space. More precisely, 
the text, we solve these equations by restricting B* ,N (p) 
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we expand B V ' N (p) in terms of orthogonal polynomial 
functions W7^(p) [see Eq. (flTJj) ] and approximate the 
series by the first three terms [see Eqs. (TTT))) and (|2"U1) ]. 
The expansion coefficients (bj^ , n = 0, 1, 2) in Eqs. (fH?)) 
and ([2H)l are determined as follows: we multiply Eq. (|B7[) 
by " ("» = M,2), E q- (ESI by 



j2^(PiPj-^r) W {L)(p) ( m = °> 1j 2 )) and integrate them 
over p. Then, by using the orthogonal conditions (|17p . 
(|18j) . we will obtain six independent equations for the 
coefficients i>7\ (n = 0,1,2). The resulting equations 
are still complicated but can be solved in a numerical 
way. 
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